Squeezed-light source for the superresolving microscopy 
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We propose a source of multimode squeezed light that can be used for the superresolving mi- 
croscopy beyond the standard quantum limit. This source is an optical parametric amplifier with 
a properly chosen diaphragm on its output and a Fourier lens. We demonstrate that such an ar- 
rangement produces squeezed prolate spheroidal waves which are the eigen modes of the optical 
imaging scheme used in microscopy. The degree of squeezing and the number of spatial modes in 
illuminating light, necessary for the effective object field reconstruction, are evaluated. 
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| Spatial behavior of nonclassical light is presently attracting an increasing interest both in theory and experiment 
Q, 0] . Quantum effects in optical imaging and other transverse spatial phenomena are studied within a framework 
of a European project "Quantum Imaging" One of the problems recently addressed in this context is about the 
■ ultimate quantum limit of the optical resolution . A classical criterion of resolution was formulated at the end of 
\^ | the last century by Abbe and Raylcigh. It states that optical resolution is limited by diffraction on the system pupil. 

However, though it might be difficult to accept, the diffraction limit is not a fundamental limit like, for example the 
' Heisenberg uncertainty relation. Recent publications have theoretically demonstrated that in principle this limit can 
^ . be beaten for both writing and reading of optical information. 

Sub-diffraction- limited optical recording of information, coined as "quantum lithography", uses a specially designed 
nonclassical light, TV-photon entangled state, together with a TV-photon absorber which allows to reach theoretically 
^-H ' a resolution of X/N @,EJ. To go beyond the diffraction limit in read-out of optical information like in imaging, is in 
, principle possible having some a priori information about the object using the so-called "superresolution" techniques. 
As shown recently in Ref. Q the ultimate limit in such superresolution is set not by diffraction but by the quantum 
fluctuations of light illuminating the object of finite size and by the vacuum fluctuations outside the object. These 
' quantum fluctuations of the illuminating light together with the vacuum fluctuations outside the object set up the 
standard quantum limit of superresolution which can be much smaller than the diffraction limit. Moreover, one can 
go beyond this standard quantum limit using a specially designed nonclassical light for illumination of the object. In 
recent experiment |7j it was demonstrated that using multimode squeezed light allows to measure a displacement of 
2.9 A of a laser beam with wavelength A = 1064 nm. 

In Ref. Kolobov and Fabre have suggested a scheme that allows to improve the superresolution beyond the 
standard quantum limit in reconstruction of an optical object using a multimode squeezed light. This scheme was 
formulated in terms of the so-called prolate spheroidal waves which are the eigenfunctions of the optical imaging 
scheme. To achieve the superresolution with multimode squeezed light one has to prepare these prolate waves in 
squeezed state. The question remained on how to produce such squeezed prolate waves. 

In this letter we provide the answer to this question. Precisely, we demonstrate that an optical parametric amplifier 
(OPA) with a properly chosen diaphragm on its output and a Fourier lens, produces squeezed prolate spheroidal waves 
used in the superresolving microscopy. We investigate the quantum statistics of the squeezed prolate spheroidal waves 
in our scheme in dependence of the physical parameters of the OPA and of the optical configuration. We formulate 
simple estimates on the number of the "object elements" to be reconstructed in connection with the number of degrees 
of freedom in the nonclassical illuminating light. 

The scheme of optical imaging with multimode squeezed light is shown in Fig. 1. For simplicity we consider a 
one-dimensional case. The part of the scheme to the right of the object plane performs diffraction-limited imaging 
of an object of finite size X located in the object plane [4|. The first lens L\ performs the Fourier transform of the 
object field into the pupil plane with a pupil of size d. After the second Fourier transform by the lens Li an image is 
created in the image plane. This image is a diffraction-limited copy of the object due to the finite size of the pupil. 

The part to the left of the object plane is an illumination scheme. It consists of a traveling-wave OPA placed in the 
source plane and a Fourier lens L. It is well-known from the literature that a traveling- wave OPA with plane- wave 
pump and nonlinear crystal with large transverse area creates multimode squeezed vacuum on its output [l|. A new 
feature of our scheme is a diaphragm of size d s on the output of the OPA which serves for selection of the transverse 
modes in squeezed state. As we demonstrate below, when the size of this diaphragm matches the size of the pupil, 
d s > d, this setup squeezes exactly the prolate spheroidal waves which are the eigen modes of the imaging scheme. 
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FIG. 1: Schematic of the optical imaging with squeezed light. 



This result can be easily understood qualitatively. Indeed, when all three lenses in the scheme have the same focal 
distance /, as in Fig. 1, the lenses L and L\ create a geometrical image of the diaphragm d s in the pupil plane. 
Therefore, it is intuitively clear that one has to match the diaphragm size d s and the pupil size d to select the modes 
of the source which will pass through the imaging scheme. 

Let us introduce the dimensionless coordinates in the object and the image plane as s = 2x/X, and the dimensionless 
coordinates in the source and the pupil plane as £ = 2y/d (see Fig. 1). The dimensionless photon annihilation operators 
in the source, object, pupil, and image planes are denoted respectively as c(£), a(s), /(£), and e(s). These operators 
obey the standard commutation relations, 

m),cHa = m-e), [aoo,aV)]=*(«-*')> w 

and similar for /(£), and e(s). Naturally, the annihilation and creation operators at different planes commute with 
each other. These operators are normalized so that {c' , for example, gives the mean photon number per unit 
dimensionless length in the source plane. Similar normalization is valid for the other planes. The Fourier transform 
(Tc) (s) physically performed by the lens L in Fig. 1 reads as follows, 

r~r~ r°° 

a(s) = (Tc) (s) = yj— J dte-^c(0, (2) 
where c = — — is the space-bandwidth product of the imaging system. 

' n 

The quantum theory of the imaging scheme in Fig . 1 was formulated in Ref. |4| in terms of the prolate spheroidal 
functions tpk(s) 0,0 (fo r some examples see Ref. .10]). These are the eigen functions of the imaging operator of the 
scheme, orthonormal on the interval — oo < s < oo. The photon annihilation operator a(s) in the object plane can be 
written as decomposition over ipk(s), 

oo 

a(s) =Y / A kMs)+N(s), (3) 

k=0 

where the operator-valued coefficients Ak are evaluated as 

/oo 
dsa(s)ij k (s). (4) 
-oo 

The operators Ak and A\ satisfy the standard commutation relations of the photon annihilation and creation operators 
for discrete modes, 

[A k ,Al,]=5 ktk ,, [A k ,A k ,} = 0. (5) 

The Fourier transform ip k (£,) of the prolate spheroidal functions tp k (s), performed by the lens Li, is zero outside the 
interval |£| < 1. Consequently, the set of functions {tpk(s)} is not complete in the Hilbert space L 2 (— 00,00), and to 
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satisfy the commutation relations £[J one has to add an additional term N(s). This term has zero Fourier spectrum 
in the interval |£| < 1 and does not contribute to the coefficients A k , 

ds^ k (s)N(s)=0. (6) 

> 

The field N(s) can be decomposed over the complementary set of prolate functions {#fc(s)}, orthogonal to {^(s)}. 

Physically speaking, the first term in Eq. @ is the object field component that propagates in our scheme through 
the pupil to the image plane. The second object field component N(s) is absorbed outside the pupil and is not 
observed. Therefore, in what follows we shall omit N(s) in the object field JHJ. 

To obtain the canonical transformation of the photon annihilation and creation operators in the imaging scheme 
on Fig. 1, one has to split the coordinates s and £ into two regions, the "core", \s\ < 1 and |£| < 1, corresponding to 
the central regions of the object (image) and the source (pupil), and the "wings", |s| > 1 and |£| > 1, outside these 
areas. The orthonormal bases in these regions of the object (image) plane are given by 



where are the eigenvalues of the corresponding prolate spheroidal functions ipk{s), depending on the space- 
bandwidth product c. We should note that the functions (fik(s) are complete in the Hilbert space L 2 (—l, 1). 
It follows from Q that the functions ipk(s) and 9 k (s) can be written in the form 



-0fc(s) = \f\kSPk{s) + V 1 - ^kXk(s), k {s) = y/l - \ k tpk(s) - y/\kXk(s)- (8) 

Similar relations take place in the source (pupil) plane. 

In terms of two sets {fk(s)} and {xfc( s )} we can write the annihilation operators in the source and the object 
planes as 

oo oo 

e(o=Ew fe (o+E4»(o. ( 9 ) 

fe=0 fe=0 

and 

oo oo 

a(s) = ^2a k ip k (s) +^2"b k Xk(s)- (10) 

fc=o fc=o 

Here c k and a k are the annihilation operators of the prolate modes (p^ in the core region of the source and the object 
planes, while dk and are the annihilation operators of the prolate modes Xk i n the wings regions. The operators c k 
and dk are expressed through the field operator c(£) by 



c k = <%c(0v k (0, d k = dtc(0xk(0- (11) 



Similar relations hold for dk, bk and a(s). 

In our analysis we shall use the following property of prolate spheroidal functions 



dSM$)e- m = (^T^^Ml/c). (12) 

Using J7J, iJSJ, the field transform J5J) between the source and the object plane and this equation, one can find the 
following propagation relations for the core and wings of the light wave emitted by the source: 



(T<p k ) (s) = {-i) k VXkipk(s) + Vi-^kXk(s) = H)Vfc(s), 



{Txh) (*) = Vl-Afe^(s)- y/X k Xk(s) =(-i) k e k {s). (13) 
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Taking into account © and l|10|) we obtain, 

oo 

a(s) = (Tc) (s) = J2(~ 

k=0 

and 

A k = {-i) k c k 



CfeV'fe(s) + d k 8 k (s) 



(14) 



(15) 



As expected, there is no contribution into A k from the second sum in Eq. © containing operators d k and describing 
the illumination coming from the wings area of the source. If the source diaphragm is larger or matches the size of 
the pupil, d s > d, it has no effect on the operator amplitudes c k and A k , see l|ll|) . and we shall neglect the diaphragm 
in the calculation of A k ■ 

To obtain explicitly the coefficients A k in case of illumination of the scheme by a traveling-wave OPA we shall use 
simplified description of an OPA with a plane- wave undepleted pump (see Ref. |l|). In this approximation one can 
find analytically the spatial Fourier amplitudes c(q), 



(16) 



of the field at the output of the crystal as a linear transformation of the corresponding input Fourier amplitudes Ci n (q) 
and c\ n (q), 



c(q) = U( q )c in (q) + V(q)cl(-q). 



(17) 



Here the operators Cm(q) at the input of the OPA are in the vacuum state. The complex coefficients U(q) and V(q) 
depend on the nonlinear susceptibility of the crystal, its length, and the matching conditions in the OPA. These 
coefficients have the property, 



\U(q)\ 2 -\V(q)f 



1, 



(18) 



that guarantees thepreservation of the commutation relations Using the standard parameters of multimodc 
squeezing (see Ref. |l|) we obtain: 



U(q)±V*(-q) 



e ±r(9) cos 6(q) +ie Tr(q) sm0(q) 



(19) 



The phase of the amplified (stretched) quadrature amplitude of the OPA output field c(q) is 0(q). The phase factor 
exp(— itp(q)) specifies the phase of the quadrature amplitudes of the input field Cj„(g) which are squeezed or stretched. 
For the vacuum input field of the OPA the last phase is irrelevant. 

The operator amplitudes A k are found explicitly with the use of Eqs. i|T5|) and lfTT|) . The quantities £„(£) and 
<fik(£,) are expressed through their Fourier transforms l|17l) and 11311 . After some calculation we obtain, 



1 



/2ttc 



dqMl/c) U(q)c m (q) + V(q)c\ n (-q) 



We shall introduce the real quadrature components of the field amplitudes A k in the object plane as 

A k = Ai k + iA 2k - 
For the variances of these quadrature components we obtain 

dqi>l(q/c) [ e ±2r(9) cos 2 %) + e^ 2r(q) sin 2 %)1 , 



((Ail*) ) 



1 

4^ 



(20) 



(21) 



(22) 



( ( Ai 2fe ) ) = ± J dq^ k (q/c) [e^M cos 2 9{q) + e ±2r ^ sin 2 9(q) 



(23) 



Here the upper and the lower sign correspond respectively to the even and the odd prolate spheroidal functions ipk- 
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It follows from this result that the prolate spheroidal waves in the object illumination can be prepared in squeezed 
state. By proper choice of the squeezing phase 9{q) at low spatial frequencies q one can minimize quantum fluctuations 
in one of the quadrature amplitudes A a k, a — 1,2 (namely, in the one detected in the image plane of our scheme). 
Taking the degree and the phase of squeezing as constant, r{q) — r, 0{q) = 9, we can estimate the variance of the 
squeezed quadrature amplitude: 

((AA ak ) 2 } ~ e- 2 74. (24) 

In reality the spatial-frequency band of multimode squeezing is limited by the phase-matching condition in the OPA. 
One has to take into account the spatial-frequency dispersion of squeezing, that is, the frequency dependence of the 
squeezing phase 9(q) due to diffraction in free space and inside the OPA. Both phenomena deteriorate squeezing of 
prolate spheroidal waves. 

As shown in 0, 0, , the effect of diffraction on squeezing can be almost perfectly compensated by means of 
adjustment of the lens array. If this is done, it follows from (I22|) . I|23(l . that the degree of squeezing of prolate spheroidal 
waves depends on the overlap in the object plane of two areas: (i) the area illuminated by the squeezed plane waves 
c(q), which are focused by lens L to the points s — q/c, and (ii) the area of the energy localization ~ ipk( s ) °f the 
prolate spheroidal waves. 

In analogy to other phenomena with spatially- multimode squeezed light, in the even and odd components of the 
object field the different quadrature amplitudes are squeezed; for, say, 9 = 0, the squeezed quadratures are An- and 
Aik for the odd and the even prolate waves respectively. In measurement this implies, e.g., the spatially resolved 
homodyne detection with numerical evaluation of the odd or even squeezed amplitudes (depending on the local 
oscillator phase). 

Finally, we can formulate the conditions on multimode squeezing in the object illumination in terms of the number 
N of independent degrees of freedom in the light field, propagating through the diaphragm of the OPA (we assume 
here the optimum size d s = d of the diaphragm). The properties of the OPA emission can be characterized by the 
coherence length l c of the output field c(£). The spatial- frequency range q c of effective squeezing is related to the 
coherence length by the estimate 

ird . . 

kd < (25) 

The minimum requirement on the OPA is that the effectively squeezed waves c(q c ) illuminate, after passing the lens 
L, the object region |s| < 1. That is, the waves 

\q c \/c>l, (26) 
should be squeezed. This gives the following estimate for the coherence length, 

«. < g, <*> 
and for the number N of independent degrees of freedom in the illuminating light, emitted from the region d: 

N ~ d/l c > ^ = S. (28) 

Here S is the Shannon number of our optical scheme. 

As seen from ©, the wave profiles ipk(s) can be expanded in terms of two bases, {ipk(s)} and {Xfc( s )}j representing 
the core and the wings of the illuminating field in the object plane, where Xk < 1 are the eigenvalues of the imaging 
transformation. As known from the theory of prolate spheroidal functions [alj, these eigenvalues are close to 1 only 
for k < 5 1 , where S is the Shannon number. For higher values of the index k the field energy in prolate spheroidal 
waves V'fe(s) is concentrated in the wings, that is, outside the object area |s| < 1. Hence, the condition N ~ S provides 
an effective squeezing only of the prolate spheroidal waves with k < S. In order to minimize quantum noise of the 
higher prolate spheroidal waves (with k > S), it is necessary to use OPA with large number of effectively squeezed 
spatial modes of radiation, 

N > S, (29) 



and to illuminate by non-classical light a spot in the object plane with the size much larger than the object itself. 
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When the condition N 3> S is met, one can illuminate the object (within the area \s\ < 1) by a bright classical 
wave with a properly chosen phase. Since for <C 1 our imaging scheme is analogous to a highly non-symmetrical 
beamsplitter p|, the squeezed illumination concentrated mainly in the wings area of the object plane allows for the 
low-noise measurement of the relevant amplitudes of the image field. 

This work was supported by the Network QUANTIM (IST-2000-26019) of the European Union and by the INTAS 
under Project No. 2001-2097. 
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